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INTRODUCTION 


Experimental  observations  indicate  that  when  metallic  materials  are 
subjected  to  cyclic  torsional  loading,  a  hardening  behavior  similar  to  the 
case  of  uniaxial  loading  occurs.  In  a  previous  paper,  Wu  and  Yip* 
successfully  applied  the  ‘endochronic  theory  of  plasticity  to  describe  the 
cyclic  hardening  phenomenon  under  uniaxial  loading,  and  it  is  the  objective  of 
this  report  to  apply  this  theory  to  the  case  of  cyclic  torsional  loading  of  a 
thin-walled  cylinder  and  a  solid  cylinder  of  circular  cross-section  and 
investigate  the  residual  stress  distribution. 

BRIEF  SUMMARY  OF  IMPROVED  ENDOCHRONIC  THEORY 

The  endochronic  theory  of  plasticity  developed  by  Valanis^  is  based  on 
the  notion  of  intrinsic  time  and  the  internal  variable  theory  of  thermo¬ 
dynamics.  The  original  definition  of  intrinsic  time  has  led  to  difficulties 
in  cases  where  the  history  of  deformation  involves  unloading.  Valanis^  has 
since  introduced  a  new  concept  on  intrinsic  time  to  overcome  these 
difficulties. 

The  intrinsic  time  C  is  a  measure  of  irreversibility.  A  definition  of 
new  intrinsic  time  which  is  more  closely  representative  of  unloading  behavior 
is  given  as  follows: 


1-Wu,  H.  C.  and  Yip,  M.  C.  ,  "Endochronic  Description  of  Cyclic  Hardening 
Behavior  for  Metallic  Materials,"  ASME  J.  of  Eng.  Materials  and  Technology, 
Vol.  103,  July  1981,  pp.  212-217. 

^Valanis,  K.  C. ,  "A  Theory  of  Viscoplasticity  Without  a  Yield  Surface,  Part  I 
and  Part  II,  Archives  of  Mechanics,  Vol.  23,  1971,  pp.  517-551. 

^Valanis,  K.  C. ,  "Fundamental  Consequences  of  a  New  Intrinsic  Time  Measure  - 
Plasticity  as  a  Limit  of  the  Endochronic  Theory,"  Archives  of  Mechanics,  Vol. 
32,  1980,  pp.  171-191. 


1 


(1) 


d dft-j^j 

dsij 

dfiij  =  deij  -  kx  -- —  (2) 

2y0 

where  U-jj  is  a  strain  like  tensor,  e^j  and  sjj  are  deviatoric  strain  and 
stress  tensors  respectively,  is  a  positive  constant  such  that  0  <  k[  <  1 


and  y0  is  shear  modulus. 


The  general  constitutive  equation  for  shear  response  of  a  material  with 


an  elastic  hydrostatic  response  and  no  coupling  between  deviatoric  and 


hydrostatic  behavior  is: 


z  deij 

sij  =  2  /  y(z-z')  -  dz'  (3) 

J  o  dz' 


where  z  is  an  intrinsic  time  which  is  related  to  C  by  the  following  time  scale 


dC 

—  -  f(c)  (4) 

dz 


It  has  been  shown  in  Reference  3  that  f(c)  describes  isotropic  hardening  and 


is  therefore  termed  the  hardening  function. 


Define 


y(z)  =  M0  G(z)  (5) 


where  G(0)  =  1  indicates  initially  elastic  response.  Using  Laplace  transform 


technique,  (for  details  see  Reference  3),  Eq.  (3)  becomes 


z  d^ij 

sii  =  2y0  J  p(z-zf ) - dz f  (6) 

J  o  dz ' 


^Valanis,  K.  C. ,  "Fundamental  Consequences  of  a  New  Intrinsic  Time  Measure  - 
Plasticity  as  a  Limit  of  the  Endochronic  Theory,"  Archives  of  Mechanics,  Vol. 
32,  1980,  pp.  171-191. 
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in  which,  for  the  case  of  k]_  *  1, 


p(z)  **  P0  <KZ)  +  PI  (z)  (?) 

where  p^(z)  is  composed  of  a  finite  sum  of  exponential  terms.  A  general 
constitutive  equation  in  terms  of  initial  yield  stress  Sy  and  plastic  strain 
^ij  (ki  =  1)  can  be  obtained  as 

dnij  z  dnij 

sii  =  sy - +  2p0  /  pi(z-z') - dz'  (8) 

J  1  dz  o  dz' 

by  substituting  Eq.  (7)  into  Eq.  (6)  and  defining  sy  =  2p0p0.  Note  that  at  z 
=  0: 


5ij 


Sy - 

y  dz 


(9) 


z=0 


Also  from  Eq.  (2),  the  condition  =  0  gives  the  relation, 

8ij  =  2y0  eij  (1°) 

Equation  (10)  merely  attests  to  the  fact  that  while  z  =  0  the  deformation 
process  is  reversible  and  therefore  the  deviatoric  stress  response  is  elastic. 


DESCRIPTION  OF  CYCLIC  SHEAR  RESPONSE 


For  pure  shear  deformation  with  Pl(z)  represented  by  one  exponential 


term,  one  has 

dx 

dft  =  dn -  (11a) 

2yc 


and 


2y0  Pl(z)  *  2 \i\  e"az  (lib) 

where  r\  is  total  shear  strain  and  ft  is  plastic  shear  strain.  Eq.  (8)  yields 

dft  rz  .  .  N  dft 

x  =  xv  —  +  2ui  /  e-ct(  z-z ) - dz’  (12) 

7  dz  o  dz' 
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with 


—  -  f(z)  and  dc  -  Idftl  (13) 

dz 

and  Ty  is  the  shear  yield  stress.  Using  Eqs.  (12)  and  (13)  and  a  suitable 
function  f,  the  governing  equations  of  torsional  test  of  a  thin-walled  tube 
during  loading,  unloading,  and  reloading  can  be  derived.  The  initially 
elastic  unloading  (reloading)  response  and  followed  Bauschinger  effect  will  be 
governed  by  the  material  property  itself,  provided  that  the  intrinsic  time 
measure  is  correct. 

If  the  first  unloading  of  stress-strain  curve  begins  when  intrinsic  time 
measure  £  reaches  t*,  the  positive  property  of  t  requires  Eq .  (13)  to  be 

de  =  -dfl  *  (14) 


Define  and  £*  as  loading  and  unloading  measures  around  the  neighborhood  of 

+ 


C*  In  t  space.  Then  Eq .  (12)  leads  to 


dC 

T  -  Ty  — 

y  dz 


* 

z_ 


dft 


+  2p  i  j  e  a(  z  z1  ) - dz 1 

dz f 


(13) 


+ 

"  Tv  - 

y  dz 


* 

z+  ,  t  d  n 

+  2m  /  e"a(z"z  ) - dz 1 

n,  dz f 


(16) 


*+ 


where  x  and  x+  are  stress  states  nearby  t*  at  loading  and  unloading  processes 


respectively.  In  the  limit,  when  q  *  and  *  C* , 

+ 

dc' 

T"  -  T+  =  2lv - 

y  dz 


(17) 
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Therefore  the  shear  stress  is  discontinuous  with  a  discontinuity  of  magnitude 


dC 

Ax  =  2xv  — 
y  dz 


(18) 


in  the  initially  elastic  unloading  region.  If  the  consequent  loading 
reversals  take  place  at‘£**,  £***,  ...  then  there  is  an  elastic  change  of 
magnitude 


(19) 

V-**  y-***  . 

c  »  c  ’  . 

in  the  shear  stress.  During  these  elastic  responses  the  values  of  C  and  z 
remain  unchanged . 

During  the  elastic  response,  dc,  is  zero  and  the  constitutive  equation  is 
simply 


d  t 

At  =  2tv  — 
y  dz 


dT  =  2p0  dh  (20) 

Once  the  change  is  larger  than  T~  -  t+,  then  the  material  behavior  is  governed 
by  Eq.  (12). 

TORSION  OF  A  CIRCULAR  CYLINDER 


In  laboratory  experiments  on  torsion  of  solid  bars,  the  recorded  data  are 
usually  the  strain  at  the  outer  fiber  and  the  amount  of  externally  applied 
torque.  In  order  to  describe  these  experimental  results,  the  shear 
constitutive  equation  established  in  the  previous  section  and  based  on  a  thin- 
walled  tubular  specimen  should  be  applied  together  with  the  equations 
discussed  in  this  section. 

The  external  torque  for  a  solid  bar  with  a  circular  cross-section  is 


given  as: 


To  =  2tt  /  rr2  dr 

3  Jo 


(21) 
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where  t  is  current  shear  stress  state  corresponding  to  location  r,  and  ra  is 
the  radius  of  cross-section.  The  torque  can  be  approximated  by  discretizing 
the  circular  cross-section  into  a  finite  number  of  concentric  circular  rings 
and  assuming  t  to  be  constant  over  each  ring.  Thus 

n 

Ts  =  2tt  l  Ti  ri2  Ari  (22) 

i=l 

Now  Eqs.  (12)  and  (13)  along  with  Eq.  (11a)  can  be  solved  to  yield  t-£  at  each 
fiber,  if  the  value  of  (shear  strain)  at  that  fiber  is  known.  Geometrical 
considerations  show  that  radial  lines  have  to  remain  straight  after 
deformation.  Thus,  one  concludes  that 

ri 

hi  -  —  na  *  (23) 

ra 

where  na  is  the  strain  at  the  outer-most  fiber. 

Recall  that  there  is  a  yield  stress  introduced  in  Eq.  (8)  when  k\  =  1  and 
C  =  0.  Hence,  an  elastic  core  always  exists  during  deformation  whose  radius 
re  is  easily  computed  as 

Ty  ra 

re  = -  (24) 

2  M0  na 

If  the  experiment  is  strain  controlled  with  strain  at  ra  varying  between  -na 
and  4na,  and  with  na  in  the  plastic  range,  then  the  elastic  core  radius 
remains  fixed  during  all  stages  of  loading  after  the  first  load  reversal. 

Now  that  is  known,  the  value  of  shear  stress  at  each  fiber  can  be 
evaluated  and  used  in  Eq.  (22).  Note  that  at  the  points  of  load  reversal,  for 
the  fibers  in  the  plastically  deformed  region,  there  exists  an  elastic  range 
governed  by  Eq.  (19).  This  discrete  type  formulation  of  torque  has  the 
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advantage  of  being  capable  of  describing  the  transient  and  residual  stresses 
in  the  bar  as  well. 

COMPUTATION  . 

The  governing  equations  are  Eqs.  (11a),  (12),  (13), 
are  summarized  below: 

dT 

dft  =  dp  -  - 

2y0 

dft  z  .  dft 

X  =  Ty  —  +  2\i\  j  e  a(z  z  ' - dz'  (12) 

dz  o  dz? 

(13) 

(22) 

ri 

^i  -  —  na  (23) 

ra 

Since  the  analytical  solution  of  these  equations  is  quite  involved,  and 
in  particular  requires  different  treatment  and  derivation  for  different 
hardening  functions  f(z) ,  a  numerical  scheme  has  been  developed  to  solve  the 
above  equations.  Since  the  relationship  between  C,  the  independent  variable, 
and  n ,  the  controlled  variable  is  indirect,  iterative  techniques  are  an 
integral  part  of  this  program. 

In  order  to  ascertain  the  degree  of  accuracy  of  numerical  methods,  a 
hypothetical  case  was  assumed.  Namely,  the  field  equations  were  solved 
analytically  using  f(z)  =  exp  (&z),  and  the  results  were  compared  with  those 


dC 

—  =  f(z)  and  dC  =  |dft| 
dz 


Ts  =  2tt  l  Ti  ri2  Ari 
i=l 


(22),  and  (23)  which 

(11a) 
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obtained  numerically.  The  difference  between  the  two  is  so  small  that  using 
numerical  techniques  does  not  introduce  any  significant  amount  of  error. 


RESULTS  AND  DISCUSSION 

Once  the  accuracy  of  the  computer  program  was  established,  it  was  used  to 
predict  the  results  of  experiments  performed  on  annealed  AISI  4142  steel.  The 
experimental  results  of  cyclic  torsion  test  obtained  by  the  Plasticity 
Research  Laboratory  at  the  University  of  Iowa  are  presented  along  with  the 
theoretical  predictions  in  Figures  1  and  2.  This  material  does  not  show  any 
appreciable  amount  of  cyclic  hardening. 

The  most  important  factor  in  theoretical  predictions  is  the  choice  of  the 
hardening  function  f(z).  In  this  computation,  the  form 


dc 

—  =  C  -  (C-l )e-$z 
dz 


(25) 


has  been  used  because  of  its  simplicity  and  its  proven  usefulness  in  cases  of 
cyclic  loading. 

Following  accepted  procedures,  the  shear  stress-strain  curve  for  the 
material  was  obtained  from  a  thin-walled  tubular  specimen,  and  the  value  of 

the  material  constant  were  determined  from  this  data.  The  values  of  the 

constants  are:  a  =  1100,  3  =  30,000,  C  =  1.9,  yQ  =  107  psi  (68.9  x  103  MPa), 

=  3.6  x  106  psi  (24.8  x  103  MPa),  Ty  =  6,500  psi  (44.79  MPa).  As  can  be 

seen  from  Figure  1  this  set  of  constants  predict  the  experimental  results 
reasonably  well.  The  same  set  of  constants  were  then  used  to  predict  the 
results  for  a  solid  bar  test.  As  evidenced  in  Figure  2,  theoretical  and 
experimental  results  are  in  reasonable  agreement. 
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As  a  consequence  of  the  computational  process  the  distribution  of  stress 
in  the  cross-section  is  evaluated.  Such  distribution  at  different  magnitudes 
of  torque  during  the  first  loading  half-cycle  is  presented  in  Figure  3. 

Notice  that  the  outer  fiber  is  the  first  one  to  yield ;  subsequently  as  more 
torque  is  applied,  the  radius  of  the  elastic  inner  core  gets  smaller.  Also 
note  that  the  rate  of  hardening  for  each  fiber  decreases  which  is,  of  course, 
in  accordance  with  strain  hardening  phenomenon. 

Figure  4  presents  the  distribution  of  the  shear  stress  in  the  bar  at 
different  stages  of  the  first  unloading.  The  flat  portion  in  the  lower  curves 
corresponds  to  the  fibers  which  have  surpassed  this  initial  elastic  unloading 
and  started  plastic  unloading.  An  interesting  observation  can  be  made  here, 
when  the  amount  of  reversed  torque  gets  larger,  the  size  of  the  plas'tically 
deformed  region,  i.e.,  the  flat  portion,  gets  bigger  as  well,  but  for  each 
individual  fiber  the  developed  plastic  stress  does  not  increase  accordingly. 
Rather,  the  applied  torque  is  preferably  compensated  for  by  elastic  relief  of 
shear  stress  in  the  internal  fibers.  In  other  words,  whenever  possible  the 
material  prefers  to  behave  elastically  rather  than  plastically.  This  is,  of 
course,  natural  because  energy  dissipation  associated  with  elastic  behavior  is 
much  smaller  than  that  of  a  plastic  one. 

The  distribution  of  residual  stresses  at  the  end  of  first  cycle  where  the 
applied  torque  is  equal  to  zero,  is  presented  in  Figure  5. 

For  purposes  of  investigating  the  implications  of  the  model  developed 
here,  a  hypothetical  material  with  appreciable  cyclic  hardening  behavior  was 
studied.  The  shear  stress-strain  behavior  of  such  material  under  fully- 
reversed  torsional  loading  is  presented  in  Figure  6.  The  material  constants 
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for  this  material  were  determined  as:  a  =  1000,  6  =  50,  C  =  1.5,  p0  =  107 
psi  (68.9  x  103  MPa),  uj  =  4  x  106  psi  (27.56  x  103  MPa)  and  xy  =  10,000  psi 
(68.9  MPa).  A  steady  loop  is  established  after  a  few  cycles. 

Figure  7  presents  the  residual  stress  distribution  at  the  end  of  reversal 
cycles  up  to  and  including  the  steady  loop. 
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SHEAR  STRESS 
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Figure  1.  Shear  Stress-Strain  Curves  for  Thin-Walled  Specimens. 
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Figure  2#  Torque  Shear-Strain  Curves  for  Solid  Circular  Cylinder, 
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Figure  3.  Stress  Distribution  at  Different  Stages  of  First  Loading  of  Solid  Cylinder. 
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Figure  4.  Stress  Distribution  at  Different  Stages  of  First  Unloading  of  Solid  Cylinder. 
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Figure  5.  Residual  Stress  Distribution  at  Zero  Torque 
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Figure  6.  Theoretical  Shear  Stress-Strain  of  a 
Hypothetical  Material. 
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Figure  7.  Residual  Stress  Distribution  in  a  Solid  Cylinder 
(Hypothetical  Material) . 
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